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Abstract 



We consider a finite-dimensional model for the motion of microscopic organisms whose propul- 
sion exploits the action of a layer of cilia covering its surface. The model couples Newton's laws 
driving the organism, considered as a rigid body, with Stokes equations governing the surrounding 
fluid. The action of the cilia is described by a set of controlled velocity fields on the surface of 
the organism. The first contribution of the paper is the proof that such a system is generically 
controllable when the space of controlled velocity fields is at least three-dimensional. We also 
provide a complete characterization of controllable systems in the case in which the organism 
has a spherical shape. Finally, we offer a complete picture of controllable and non-controllable 
systems under the additional hypothesis that the organism and the fluid have densities of the 
same order of magnitude. 

1 Introduction 

The mathematical description of the motion of swimming organisms and of the mechanisms of their 
propulsion is a challenging and wide field of research. This paper aims at giving a contribution to 
such a domain in a comparatively narrow direction. 

The object of our study are microscopic organisms. The most prominent aspect of the models 
describing their motion is the very high level of viscosity ([3l [T8l I22j). The mechanism of their 
propulsion depends heavily on the different species and admits a great variety of models (see, for 
instance, [21 [71 [H [T71 and references therein) . 

Among microscopic organisms we will further restrict our attention to the class of ciliata, whose 
propulsion is determined by cilia, hair-like organelles covering their surface and whose size is very 
small compared with that of the organism (see [11 [H [9]). Recently San Martin, Takahashi and 
Tucsnak [19] proposed a control theory approach to the description of ciliata. Following the former 
literature (and in particular [9l [16]) they assume that the organism is a rigid body and that the 
motion of the cilia is described by a set of controlled velocity fields on a surface enclosing the layer 
of cilia. They are therefore faced to an infinite-dimensional control system coupling Navier-Stokes 
equations, which describe the motion of the fluid, and the equations for the rigid body arising form 
Newton's laws. The hypothesis of high viscosity (ie, of low Raynolds number) mentioned above is 
exploited to reduce such a control system to a finite-dimensional version of it. The state-space is 
given by the coordinates of the center of mass of the organism, its orientation and its angular and 
linear velocities. 

Summarizing, San Martin, Takahashi and Tucsnak single out a 12-dimensional nonlinear control 
system which is affine in the controls. The control system depends on a number of parameters which 
are explicitly derived as functions of the shape of the organism, of its mass distribution, and of the 
surface velocity fields describing the motion of the cilia. The authors study the control properties 
of the system proving that it is generically controllable when the number of controlled vector fields 
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is at least six (see Section [5] for details). They also provide some example of controllable systems 
in the case in which the organism has a spherical shape. In this special case they also study the 
controllability properties of the linearization of the system at its rest position. 

Our scope is to enhance the results of [19] basically in three directions: firstly we want to give 
sharper results on the generic controllability of the system; secondly we want to obtain a complete 
characterization of controllable systems in the case in which the organism is spherical; thirdly we want 
to provide a complete picture of controllable and non-controllable systems in the simplified situation 
which -following again [l9]- comes into effect in the case in which the organism and the fluid have 
densities of the same order of magnitude. The overall aim is to provide a better understanding of the 
dependence of the controllability properties of the model proposed in [19] on the physical quantities 
characterizing it. The hard aspect of the model is indeed the difficulty to obtain reliable expressions 
for the surface velocity fields describing the motions of the cilia. In this regard, it is important to 
provide as precise as possible parametric analysis of the controllability properties of the system. 

The paper is organized as follows. Section [2] introduces the control theoretical language adopted 
throughout the paper. Section [3] presents the control system obtained in [19] and recalls the rela- 
tions between the physical quantities characterizing the organism and the parameters appearing in 
the model. Section U] is technical and contains an abstract controllability result used in the later 
sections. Section [5] is devoted to the analysis of the genericity of the controllability of the system 
introduced in Section [3j Genericity is initially formulated in terms of the parameters appearing 
in the finite-dimensional system and then expressed in terms of the physical objects they depend 
on (Section 15. 3p . The main result is that the system is generically controllable if the number of 
controlled vector fields is at least three and that one controlled vector field suffices for its generic 
accessibility (Theorem 15. ip . Section [6] analyzes the case in which the organism is spherical and 
provides a complete characterization of the sets of parameters that make the system controllable 
(Theorem 16.10 . Finally, Section [7] establishes a characterization of controllable systems when the 
organism can be assumed to have the same density as the fluid (Theorem 17. 2p . The result follows 
from a further reduction of the system that transforms it into a six-dimensional one. The reduction 
procedure, suggested in [19], is carried out in full details. 



2 Notations and definitions 

Let us first introduce some notations: in the following Id denotes the 3x3 identity matrix, Oyyixn 
is the zero m x n matrix, 0^ = Omxi is the zero vector in R"^. We denote by ^mxn the set of 
m X n real matrices. The interior of a set W is denoted by Int(14/') and its closure by Clos(M^). For 
simplicity of notation, when no confusion is possible, row and column vectors are identified. The 
canonical orthonormal basis of R'^ is denoted by {61,62,63}. For every uj £ R^ we denote by S{uj) 
the skew symmetric matrix such that ^(ci;)^ = uj x z for every z € R^, that is, 

/ -tus L02 \ 
S{uj) = { UJ3 -uji . 

\ — UJ2 / 

In what follows we deal with control systems of the type 

q = f{q,u), qeM, n G R"^, (1) 

where M is a smooth (C°°) manifold and f : M x R"^ — > TM is smooth with respect to both variables. 
Moreover, for every T > the endpoint mapping Et which associates to an initial condition go G M 
and a control function u S L°°([0, T], R™) the final point of the corresponding trajectory of ([T|) is 
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well defined and continuous on M x L°°{[0, T], R"*), where the norm is chosen on L°°{[0, T], R™). 
(For these and finer properties of the endpoint mapping see, for instance, [1].) 

Let ^ = {/(■,«) I u € R™} be the family of vector fields characterizing (pQ). For every X = 
/(•, n) S ^ we denote by e^-^ the one-parameter group generated by X, that is, e^^{q) = Et{q, til[o,j]). 
The attainable set from (7 G M is the set 

A{q,^) = UT>oET{q,L^{[0,T],Kn). 

When no confusion is possible we write A{q) = A{q, 

We say that ([T]) (or, equivalently, ^) is controllable if for every q £ M the set A{q) is equal to 
M. We say that ([T]) is approximatively controllable if A{q) is dense in M. 

Give two vector fields X and Y on M, the Lie bracket between X and Y is defined, in a local 
system of coordinates, by the relation 

[X, Y] = (DY) X - (DX) Y 

where DX denotes the derivative of X. Denote by L the Lie algebra generated by System ([1]) is 
Lie-bracket generating if for every q £ M the set L(q) = {V{q) | ^ G L} is equal to TgM. 



3 The dynamics of swimming microscopic organisms 

In [19] the authors describe the swimming of a microscopic organism immersed in an infinite volume 
of fluid by coupling the dynamics of a rigid body (representing the organism) with the Navier-Stokes 
equations describing the behavior of the fluid outside the body. Exploiting the hypothesis of very 
low Reynolds number, they derive a simplified version of the coupled system. The system obtained 
in this way is finite-dimensional and nonlinear; more precisely, it has the following expression: 

z = Az + E{z) + Bu, (2) 

C = RC, (3) 

R = RSiio), (4) 
where z = (^,w) G R^ x R^, ( G R•^ R ^ SO{3), and 

^ ' \ J [{Juj) ^ UJ) J 

Here J G ^3x3 denotes the inertia matrix of the rigid body representing the organism. Recall that 
J is symmetric and positive definite. The matrix A G ^6x6 is a function of the rigid body and the 
viscosity of the fluid and its expression is given at the end of the section. 

The control function u takes values in R"*, m > 1, and is a 6-by-m matrix, which we write as 




where both Bi and B2 belong to ^sxm- 

An important feature of the system above is that ([2]) is a well-defined control system on R^ , since 
it does not depend on nor R. 

Denote by Xq the vector field on R^ x R'^ x S0(3) which is the drift of the control system ©-([J]), 
that is, 

/ Az + E{z) \ 
Xo{z,C,R)=\ R^ 

\ RSiuj) J 
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Let moreover Xi, . . . , be the controlled vector fields, so that system ©-([H) can be written as 

m 

q = Xo{q) + ^^X^{<l), g € X 50(3). 

i=l 

Define L as the Lie algebra generated by Xq, Xi, . . . , Xm- 

We say that a vector field X on x S0{3) is a constant vector field if there exist vi,V2 G R"^ such 
that X : i-^ ivi,V2,06) and we write X = {vi,V2,0q). Clearly, Xi, . . . , Xm are constant 

vector fields. Notice that ii V = {vi,V2,0q) and W = {wi,W2,0q) are constant vector fields, then 

(V2 X Wi + W2 X Vi \ 
J-l {JW2 XV2 + JV2 X W2) (5) 
06 J 

is a constant vector field as well. 

Remark 3.1. For every constant vector field V = (^"1,^^2)06), the bracket [Xo,V] is of the form 

[Xo, V]{z, C, R) = {G{z, V), -Rvi, -RS{V2)) 

for some smooth map G with values in R^ . Therefore, if the constant vector fields contained in L 
span R^ X {Og}, then system (^-^ is Lie-bracket generating. 

We conclude the section by recalling the analytic characterization of the coefficients appearing 
in A, B and J. Let Q C R^ be an open bounded connected set whose boundary is of class C^, 
representing the (rigid) shape of the organism. Then the coefficients Jij of J are given by 

Jij = / S{x){ei X {x - ^)) • {cj X {x - 0)dx, (6) 
Jn 

where ^ denotes the mass center of Q and 6 : 0, ^ {0, oo) is the (possibly non-constant) density 
function. 

Denote by F the set R^ \ Clos($7) and, for / G N, 1 < g < oo, let 

Z)''9(F) = ye Ll^ I d'^f e L'i{F) for all a G N^, |a| = /}. 

According to |19l Lemma 3.1], for every i G {1,2,3} there exists a unique solution {h^^\p^^^) of 

-A/i(*) + Vp^^ = on F 
div/i(^) = on F 



= a on dn 
lim|^l_^/iW(y) = 



and a unique solution (ij(*),p(*)) of 

_Ai7(^) + VP(^) = on F 

divH^^ = on F 

H^^\y) = Gi X y on 

lim|^|^^i/»(y) = 

such that (£ L%F) n D^''XF) n D'^^^{F) nC^^iF) andj9«,pW G L'"(F) n Z)i'^(F) n C°°(F) 

for s G (3, oo], r G (3/2, oo] and 9 G (1, oo). 
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The Cauchy stress is the tensor a defined by the relation 



\dyi dykj^<i^k<3 



where // is the Reynolds number of the fluid. For i = 1,2, 3, define 
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where n is the unit inner normal to dQ. 

Denote by G^, 0^, T^, the matrices in ^3x3 defined by 



dQ 



Th = - / g'^Us, 



gfds, el = - (xx g^'A ds, 
[x X G^*)) ds, 



'an Jan ^ 



where is the surface element of dfl. It turns out (see [TT]) that = {Q'^)^ and that 

where m denotes the mass of the organism. Moreover, with respect to the inner product in 
defined by 

3 3 

(a, b)j = m^aibi + ^ Jika^+ibs+k, (9) 
1=1 i,k=i 

A is self-adjoint and negative-definite. 
The matrix B is defined by 

where the entries of k and K, are 



tii = - / 9^^ ■ i^jds, ICij = - G^^ ■ il^jds, 
Jan Jan 

and ipi, . . . , tpm are fixed functions in C'^{dQ, R'^), each one associated to a component of the control 
u, that describe the admissible propulsive actions of the organism. 

4 A controllability result in a more general setting 

The following proposition provides sufficient (and necessary) conditions for the controllability of 
systems that generalize those introduced in the previous section, in the sense that a more general 
structure for the dynamics in the coordinates z is allowed. 

Proposition 4.1. Consider the following control system 

z = f{z,u), zGR^ uGR"". (11) 

Then [^).^). [Tl\) is controllable if and only if it is Lie-bracket generating at {z,(,R) = (09,Id) and 
system is controllable. 
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Proof. One direction of the equivalence being obvious, let us assume that system (jlip is controllable 
and that ([3]).(j^. pT]) is Lie-bracket generating at {zX,R) = (09,Id). 
Notice that 

{z',C',R') ^ Az,C,R) =^ {z',( + RC',RR') eA{z,C + RC,RR) (12) 

for every ((^, R) G x SO {3), where the letter A is used here to denote attainable sets for system 
([3]).(|^. pT]) . Indeed, if t u{t) is an admissible control steering {zX,R) to {z'X',R') and t i— > 
{z{t), ({t), R{t)) is the corresponding trajectory, then t ^ {z{t),( + R({t), RR{t)) is an admissible 
trajectory corresponding to the same control u. 

Fix {zoXoi Ro)j {zi,CijRi) £ R-^ ^ S0{3). Since (fTTI) is controllable, then there exist ryo,?7i G 
and Qo,Qi G 50(3) such that (06,r7o,Qo) S ^^(zo, Co, -Ro) and (2;i,Ci,-Ri) G ^(Oe, r/i, Qi). We are 
left to prove that -4(06, Qo) contains {Og} x x 50(3). 

Since ([3|) . (j^ . (fTT|) is Lie-bracket generating at {z,C,R) = (09,Id), then the set 

G = Int(^(09,Id)) 

is nonempty, as it follows from Krener's theorem (see, for instance, [12} Theorem 1, p. 66]). Fix 
(2:*, C*, R*) £ G. Since (fTTI) is controllable, there exists a control law t 1— > u{t), defined on an interval 
[0,T] and such that the solution of (jlip with initial condition z{0) = z^, satisfies z{T) = Oq. The 
flow corresponding to the control law u, evaluated at time T, is a diffeomorphism sending (z*, C*, R*) 
to a point of the form (Og, C**) -R**)- particular, since the image of G by such diffeomorphism 
is contained in G, there exists an open nonempty subset V of R'^ x SO (3) such that {Og} x 1/ is 
contained in G. 

Notice the following consequence of (fT2]) : if {Oq,(jR) £ -^(Oe, ??0i Qo) and {Og,C\R') ^ C, then 

(Oe, C + ^C, RR') G Int(^(06, r/o, Qo)) • (13) 

Indeed, taking {z,C,R) = (09,Id) in (112p . we have that for every {z" ,C," ,R") in a neighborhood of 
(06,C',^')> the point {z" X + RC" , RR") belongs to A{'dQ,C,,R) and, since {^&,CR) G ^(Oe, Qo), 
we are done. 

Notice that for every nonempty open subset O of 50(3) 

W = {P1P2 • • -Pfc I A: G N,Pi, . . . ,Pfc G O} (14) 

is equal to 50(3). Indeed, since O contains at least one element Pq of finite order, ie, an axial 
rotation of angle commensurable with vr, then Id belongs to the interior of W . The completeness of 
50 (3) guarantees that W = 50(3). 

Take as O the projection of V on 50(3). Then for every R G 50(3) there exists C, G R^ 
such that (06,C,-R) belongs to O; this can be seen by noticing that there exist k elements 
(Oe, Ci, -Pi), • • • , (Oe, Cfc, -Pfc) in V such that Pi £ O and P1P2 ■ ■ ■ Pk = R and by applying repeat- 
edly (fT3]) in the special case 770 = O3 and Qq = Id. 

As a consequence, without loss of generality, 

Qo = Id, Qi = Id. 

Indeed, since there exists Co S R-^ such that (06,Co,Qo^) S C, equation (fT3l) with (C, -R) = (^o,Qo) 
shows that (Oe,^ + QoCo,Id) G ^(06,'?/o,Qo); similarly, the existence of Ci £ R-^ such that 
(06,Ci,Qi) € G implies that (06,r/i,Qi) G ^(06,r/i -Ci,Id). 

Moreover, for every element p of RP^, the Grassmannian of one-dimensional subspaces of R'^, 
there exists Cp G H-"^ such that (Og, Cp, -Rp) lies in G, where iip denotes the rotation of angle vr around 
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the axis p. Applying (fT3|) to (C;-R) = (C'l^O = {Cp^Rp) ™ the special case r/o = O3, we obtain that 
(06,Cp + -RpCpjId) lies in G. Notice that Cp + RpCp= '^{Qp^p)'P belongs to the axis p. Therefore, G 
intersects {Og} x p x {Id} for every p € RP^. 

Of special interest will be for us the expression of (jl3p when R = R' = Id, namely 

(06,C,Id)G^(06,r?o,Id), (06,C',Id)eG =^ (Og, C + C', Id) G Int(^(06, Id)) . (15) 

Let 

Sl = {v(^S'^\ {Og} X R>o(±t;) X {Id} n G / 0} 

where 5^ denotes the unit sphere in R^. As we remarked above, S\ and 5^ cover 5^ and, by 
construction, they are open. Moreover, both are nonempty, since v € implies that —v S 5^. 
Since 5^ is connected we have 5^ H 5^ ^ 0, ie, there exist v € R^ and two non-negative constants 
Ai,A2 such that both (06,Aif,Id) and (Oe, — A2'y, Id) lie in G; moreover, due to the openness of G, 
these constants can be assumed to be positive and commensurable and so there exist two positive 
integers ni and 77,2 such that niAi — n2A2 = 0. Applying repeatedly (fT5]l in the special case % = O3, 
we obtain that G contains (09,Id) and, therefore, {Og} x R^ x {Id}. 

Applying (fT5]) again with ^ = % we deduce that (06,?7i,Id) is attainable from {%,rjQ,ld). □ 

5 Generic properties 

The expression generic is commonly used to denote a property of a system that is, in a suitable 
sense, stable under small perturbations and that, even when it fails to apply, can be forced to hold 
by applying to the system an arbitrary small perturbation. 

In order to define precisely what a generic property means in the present context, define 

rJ^' = {{A, B, J) e ^6x6 X ^6xm X ^3x3 \ J = > 0, Ag Symj}, 

where 

SymJ = {A £ ^6x6 I A is symmetric and negative definite with respect to (•, •) j} 

and (•, •) J is defined as in ([9]). 

We say that a property is generic for system when m = rriQ if there exists an open and 

dense set O in H^J^" such that the property holds for every system defined by a triple in O. 

Such a definition of genericity, given in terms of the entries of the matrices A, J and B, is adapted 
to the finite-dimensional formulation ([2|)-([4|) but can appear unsatisfactory from the point of view 
of the physical model. Section 15.31 discusses how to define the genericity directly in the space of 
configurations of the organism and provides a physically justified counterpart of the following result. 

Theorem 5.1. (a) Generically when m > 1 the rank of L is maximal at every point, (h) Generically 
when m > 3 system (f^)-® ^-^ controllable. 

The generic conditions ensuring the Lie-bracket generating condition and the controllability of 
©-(II]) are given explicitly in the Propositions 15.21 and 15.31 Notice that the techniques applied in 
this section are independent of the structure imposed on A by the physical motivations of the model, 
namely, the fact that A G SymJ". Indeed, the conditions thatare obtained in the following sections 
define an open and dense set in 

E'P = {{A, B, J) G ^6x6 X ^6xm X ^3x3 \ J = > 0} 

for m > 1 (for the Lie-bracket generating condition) or m > 3 (for the controllability) whose in- 
tersection with is open and dense in H™. In other words. Lie-bracket generating condition and 
controllability are generic properties (for m = 1 and m = 3, respectively) also in the class of systems 
of the form ©-([I]) for which the triple {A, B, J) is taken in S™. 
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5.1 Generic Lie-bracket generating condition for m = 1 

When m = 1 the two matrices Bi and B2 are column vectors that we denote by bi and 62 respectively. 
Then Xi = (61, 62, Og) and L is the Lie algebra generated by Xq and Xi. 

Proposition 5.2. Let i = or i = 1. There exists an open and dense set of triples (A, B, J) in H| 
such that system (^-^ is Lie-bracket generating. More precisely, if J, bi, and 62 satisfy the three 
following conditions 

det{bi,Jb2,J%) + 0, (16) 
det(62, ^^62) + 0, (17) 

then L is maximal at every point. 

Proof. Let us prove, first of all, that if C G ^(3 x 3) is symmetric and invertible, and if 6 G is 
such that Cb and C^b are linearly independent, then 

C-\Cbxb) =det{C-^){CHxCb). (19) 

Notice that {C~\Cb x b),Cb) = = {C~\Cb x b),C^b), and therefore 

C-^{Cbxb) = k{CHxCb), (20) 

for some /c G R. For every x G R^ we have det(C~'^x, Cb, b) = (det C~^) det(2;, C^6, Cb) and also 

dei{C-^x,Cb,b) = {C-^x,Cb x b) = {x,C-\Cb x b)) = k{x,CH x Cb) = k det{x, C\Cb). 

Therefore 

(det C~^) det(x, C^b, Cb) = k det{x, C^b, Cb) 
for every x G R^, proving that (jl9p holds true. In particular, taking J = C and b = b2, we have 

J~^{Jb2 X 62) = det(J-^)(j262 X J&2). (21) 

Denote by Vi the Lie bracket [[Xq, Xi], Xi]. According to ^ we have 

Vi = 2{{b2xbi),J'\jb2xb2),0e). 

Similarly, the definition V2 = -[[Xq,Vi],Vi] leads to the expression 

V2 = (- {j~Hjb2 X b2),b2) bi + {J~\jb2 X b2),bi) b2,-{j~\jb2 X 62), 62) 62,05). 

Since {J~^{Jb2 x 62), 62) Xi + V2 = {{J~^{Jb2 x 62), 61) 62,09) we have that 

Zi = (62,09) 



belongs to L, as it follows from (jl6p and from (j2ip . 

Let Z2 = 2 [["'^O' "^i]' ^i]" "'^PP^yi^S again ([5]), we have that 

Z2 = (-62 X J-i(J62 X 62), O9). 
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Equation ([2T|) implies that 

Z2 = ((det J^i)62 X (J62 X J%),09). 

Letting 

63 = 62 X (J52 X 7^62) 
we can represent Z3 = 2[["'^0' as 

Z3 = (62 X 63,09). 
We have shown that the Lie algebra L contains the vector fields 

Zi = (62,09), Z2 = (63,09), Z3 = (62 X 63,09). 

Notice that 

63 = (62, J%) J62 - (62, J62) 7^62 

is not in span(62), since (62, J62) > and because of (fT7|) . As a consequence, the vectors 62, 63 and 
62 X 63 are linearly independent. 

Therefore, L contains every constant vector field of the type (^,09), with v in R3. It follows that 

Wi = (03,62,06) 

is in L. As a consequence of ([5]), [[Xq, W^i], VFi] = 2(03, J~^(J62 x 62), Oe) lies in L as well. Due to 
(pT|) . we deduce that 

W2 = {03,J^b2 X J62,06) 

belongs to L. Let 

64 = J^62 X J62 

and let W3 = [[Xq, M^i], W2]; we have W3 = (03,65,06) with 65 = J'^iJb^ x 62 + J62 x 64). We 
are left to prove that the vectors 62, 64 and 65 are linearly independent. Let D = det(62, 64, 65) and 
notice that 

D = det(62, 64, J"^(J64 X 62 + J62 x 64)) 

= (det J~^) det( J62, J64, J64 x 62 + J62 x 64) 

= (det J"^) ( J62, J64 X ( J64 X 62) + J64 X ( J62 X 64)) 

= (det J~^) ( J62, (J64, 64) J62 - (J64, J62) 64 + (J64, 62) J64 - II J64f 62) . 

Since ( J64, J62) = (64, 7^62) = and ( J64, 62) = (64, J62) = we have 

D = (det J-i) ( J62, (J64, 64) J62 - II J64f 62) 
= (det J-i)(||J62f (764,64) - ||J64f (J62,62)). 

Now, taking C = J'^ and 6 = 7^62 in (fT9l) . we have J64 = (det J){Jb2 x 62), from which we deduce 

(J64,64) = (detJ)((J62,62)(J62,j'62>- ||J62f). 

We also have 

||J64f = (det J)2||J62 X 62^ 

= (detJ)2(||J62f||62|p-(J62,62)'). 
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Defining A = II J62IP/ (^^2,^2) we get 

D = \\Jb2fi{Jb2,b2) {J%,Jb2) - \\Jb2\\^) - (det J) {Jb2,b2) (11^62^ 1162^ - {Jb2M?) 



,3 / \J^b2M) . .3 ,^ f \\b2f 



Cay ley-Hamilton theorem implies that 

= 02 + ai J + (det J)Id, oi, 02 £ R-i 



and thus 



[j^b2,Jb2) _ {J^b2M) , ,^ , \\b: 



^'i~rFl — TV + 01 + (det J)- 
02A + oi + (det J)- 



'21 



|2 



(J62,^2) (Jfe2,^'2) ^ '{Jb2,b2) 



{Jb2,b2 



Therefore, 



D = ( J62, 62)^ ( (a2X + ai + (det J)J^^) X - - (det J) f yJ^^A - 1 

VV {J 02,02) J \{Jb2,02) 

= {Jb2, b2f (-A^ + a2A2 + oiA + det J) 
= (J62,fe2>^det(J- Aid) 

which is different from zero because of (1181). □ 



5.2 Generic controllability in the case m = 3 

The aim of this section is to prove the genericity of the controllability of ([2]) in the case m = 3. Let 
us introduce the notation 

'Au Ai2' 
.A21 A22, 



A 



with each Aij belonging to >^3x3- 

Proposition 5.3. Let i = or i = 1. There exists an open and dense set of triples {A,B,J) in 
such that system ^ is controllable. More precisely, if B2 and An + A^^ are invertible with 
All = All — BiB2^A2i and if at least one eigenvector of J is not an eigenvector of BiB2^ , then 
system ^ is controllable. 

Propositions 15.2 1 and 15.3] lead to the following statement, which follows from Proposition 14. 1[ 

Corollary 5.4. For m > 3 and i = or i = 1, there exists an open and dense set of triples {A, B, J) 
in such that system is controllable. 

The proof of Proposition 15.31 will be split in several steps. The following lemma allows us to study 
the controllability of by investigating the controllability of an equivalent system in R'^ (instead 
of R6). 

Lemma 5.5. Assume that B2 is invertible and define B = BiB2^ and An = An — BA21. Then, 
the controllability of ([2]) is equivalent to the controllability of the system 

X = Aiix - V X X + {AiiB + A12 - BA22)v - V X Bv - BJ~'^{Jv X v), x G R^, (22) 

with the control v taking values in R"^. 
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Proof. Suppose that system ([2]) is controllable. Let x = S. — Blo and rewrite system ([2]) as 
\ X = A\\x — (J X X + {A\\B + A\2 — BA22)uJ — uj x Buj — BJ~^{Juj x uj), 



A21X + {A21B + A22)io + J~^(Ja; Xio)+B2 u. 



Given two points (xo,u;o) and in R^, there exists a control u : [0,T] — > that steers 

(xqjWo) to {xi,uji). Thus, denoting by {x{-),uj{-)) the corresponding trajectory, the control v(t) = 
uj{t) steers system (p2|) from xq to xi. 

Conversely, let system (j22p be controllable and fix two pairs (CojI^o) ^-^id in R^. Notice 

that p2p can be controlled by smooth controls u : [0, T] ^ R'^ satisfying 

v{Q)=uJo, v{T)=ui. (23) 

(The result follows from the controllability of (j22p and the density of {v € C°°([0, T], R^) | v{0) = 
cuq, v{T) = wi} in L°°(0, T) with respect to the L^-norm. The proof can be deduced from the general 
results in [10] or easily adapted from [6] and [12^ Theorem 4, p. 110].) 

Therefore, if i; G C°°([0, T], R^) satisfies ([23]) and steers ^ from ^0 - Bujq to ^1 - Bloi, then 

u(t) = 52^1 (v{t) - A2ix{t) - {A21B + ^22)^^(0 - J'\Mt) X v{t))^ 

steers ([2]) from (^Oji^o) to (^i,u;i). □ 

The following lemma establishes the controllability of system ()22p in a first case. 

Lemma 5.6. Assume that the matrix An + Aj^ has two nonzero eigenvalues with opposite signs. 
Then there exists an open and dense subset O of {{B, J) G ^6x3 x -^3x3 | J = J""" > 0} such 
that system (f22l) is controllable if {B,J) £ O. More precisely, if B2 is invertible and if at least one 
eigenvector of J is not an eigenvector of BiB2^ , then system (j22p is controllable. 



For the proof of this lemma we shall use the following result, proven in [15j . about the control- 
lability of a family of affine vector fields in R". An affine vector field X is a mapping from R" to 
R" of the form X : x Cx + c where C G ^nxn and c is a constant vector in R"; the linear part 
of X, denoted by A(X), is the linear vector field x 1— > Cx. If ^ is a family of affine vector fields, 
we denote by A(^) the set of linear parts of the vector fields in ^ and we say that ^ has no fixed 
point if there does not exist a point xq € R" such that X(xo) = for every X G 

Theorem 5.7 (Jurdjevic and Sallet). Let ^ be a family of affine vector fields given on Tl^. Assume 
that .a^ has no fixed point. If is controllable on R" \ {0} then is controllable on R". 

Proof of Lemma 15.61 Prom Lemma 15. 5|, we know that it is sufficient to prove the controllability of 
the system defined by the family of affine vector fields ^ = { C{v)x + c{v) \ v G R'^ } where 

C{v) = An - S{v), c{v) = {AiiB + A12 - BA22)v - v x Bv - BJ-^Jv x v). 

Choose fo in such a way that S{vo) = {An - Aj^)/2. Then 

An + 



111 



S{vo) 



2 

Define the family of linear vector fields ^ = { An — S{vo), S{v) \ v G R^ }. The two closed convex 
cones generated by the families A(^) and are identical, since for every finite family (ai)o<i<Af of 
positive numbers and every choice of vi, . . . ,V]\f G R^, 

N ^ N 

ao{An - S{vo)) + y^aiS{vi) = ao{An - S{vo)) + lim - qJJlh - ^(-r Vj)) 

^ — ' — >+oo r ^ — ' 

i=l i=l 
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and 



N / ^ \ ^ 

^ai{Aii - S{vi)) = {^aij [An - S{vo)) + ^aiS{vo - Vi). 

i=l \i=l / i=l 

Therefore, the controllabihty of A(^) on \ {0} is equivalent to the one of 5f (see \13\ E]). 

The trajectories of the vector field S{v) are circles contained in planes orthogonal to v and whose 
centers are at the intersections of these planes with the line Hv. Let x € R'^ \ {0}; thanks to the 
vector fields S(v), the attainable set from x for ^ contains the sphere of center O3 passing through x. 
Now, as ^11 — S{vo) has two eigenvalues with opposite sign, we can move (thanks to this vector field) 
along a direction towards the origin and along an half-line exiting the sphere and going to infinity. 
Therefore, using again the fact that the family {S{v) \ v G R'^} is transitive on every sphere, we 
proved that and thus A(^), is controllable on R^ \ {0}. 

Assume now that we can find a fixed point, denoted by xq, which is common to all the vector 
fields in For every vector v € R^, we have 

C{v)xo + c{v) = C{-v)xo + c{-v) = O3 

from which we deduce 

Aiixo-v X Bv - BJ'^{Jv xv) = 0. (24) 
Taking f = O3 gives AuXq = O3 and therefore (I24p can be rewritten as 

Bv X V = BJ-^{Jv X v) (25) 

for every v G R'^. In particular, if v is an eigenvector of J but not of B, then the right-hand side of 
(f25|) is equal to zero, while the right-hand side is not, leading to a contradiction. □ 

We are left to deal with the case where the eigenvalues of 

An + 



A 
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the symmetric part of An, are all positive or all negative. Although in this case the linear part of 
the family ^ is not controllable on R^ \ {0}, we can nevertheless adapt the method introduced in 
|15j in order to prove the controllability of (j22p . 



Lemma 5.8. Assume that the matrix Af^ is invertible and that its eigenvalues have all the same 
sign. Then there exists an open and dense subset O of {{B, J) € ^6x3 x -^3x3 I J = J"^ > 0} such 
that system (|22p is controllable if {B,J) G O. More precisely, if B2 is invertible and if at least one 
eigenvector of J is not an eigenvector of B, then system ()22p is controllable. 

Notice that, together with Lemma 15.61 Lemma 15.81 concludes the proof of Proposition 15. 3i In 
order to show Lemma 15.81 we need some definitions and have to prove some intermediate results. 

Let ^ and A(^) be defined as above. In the proof of Lemma 15.61 we pointed out that starting 
form a point x € R^ and following all possible linear vector fields of the type S{v) one can attain the 
entire sphere of center O3 and radius ||2;||. If the eigenvalues of Af^ are negative (resp. positive), the 
vector field Af^ points towards the interior (resp. the exterior) of this sphere, from which we deduce 
that 

{x} U {y e R3 j < ||y|| < ||x||} if if^ < 0, 
{x} U {y G R3 I ||y|| > ||x||} if if ^ > 0. 

In the sequel, we shall need the notion of normal accessibility, which is recalled below. 



A{x,A{^)) 
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Definition 5.9. Let Y he a family of complete vector fields on a manifold M. The point y is said 
to be normally -accessible from x if there exist , . . . , in and ti, . . . ,tp > such that 

and the mapping 

defined in a neighborhood of (ti, . . . , tp), is of rank equal to dimM at (ti, . . . , tp). 

In [15] the authors exploit the fact that if a family of vector fields is controllable, then every point 
y is normally accessible from every point x. As the family A(^) is not controllable on \ {0}, we 
are lead to prove directly the following normal accessibility property. 

Lemma 5.10. Assume that the eigenvalues ofAu+Aj-^ have all the same sign. Then we can extract 
from ^ a finite family such that there exists a sphere of center O3 and radius r of points normally 
accessible from every point of the unit sphere. 

Proof. Let vq € be such that A\^ = — S{vq) = A\^. Let, moreover, wi be an eigenvector of 
A\^ and take W2,w^ € R"^ such that the family {wi,W2,W'i) is an orthogonal basis of R^. Consider 
now the family 

A(^0) = - k2S{w2),A\^ - ksSiws) } 

extracted from A(^). 

Since S{w2) and S{w3) are nonzero skew-symmetric matrices, each of their spectra contains the 
value zero and two nonzero purely imaginary eigenvalues. Notice that 

Au-kjS{wj) 
lim — = —S{wj). 

k j — ^ -|- 00 fCj' 

As the roots of a characteristic polynomial depend continuously on the coefficients of the related 
matrix, it follows that for /c2 and k^ large enough Af^ — k2S{w2) and ^f;^ — ksS^ws) have a pair 
of non-real eigenvalues. By projecting the vector fields — k2S{w2) and A\i — k^S^w^) on the 
unit sphere, it is easy to see that the family A{^^) is transitive on the directions, that is, from 
any half-line starting from the origin, one can reach any other half-line. Moreover, thanks to the 
eigendirection wi of the vector field Af^ we can go as far as we want (case A^^ > 0) or as close as 
we want to the origin (case Al^ < 0). Hence, for every x € R^ \ {O3}, A{x, A{^^)) contains a set of 
the form R^ \ B{0,r) (case if^ > 0) or B{0,r) \ {O3} (case if^ < 0). The interior of A{x,A{^°)) 
being nonempty, there exists at least one point, denoted by y, that can be normally A(^'^)-accessed 
from X (see [21] )• Notice that since y is normally A(^'^)-accessible from x, then y is also normally 
A(^'^)-accessible from every point in a (sufficiently small) neighborhood of x. Moreover, all the 
points in A{y, A(^^)) are normally A(^'^)-accessible from x. 

We conclude by using the compactness of the unit sphere: for every x € 5^, there exists a 
neighborhood Vx of x and a set R^ \ -6(0, rx) (case A^^ > 0) or B{0, rx) \ {O3} (case Af^ < 0) whose 
points are normally A(^'^)-accessible from every point of V^. As we can include the unit sphere in 
a finite union of neighborhoods Vx, we can claim the existence of a sphere of radius r whose points 
are normally A(^^)-accessible from every point of the unit sphere. □ 

The following lemma guarantees the unboundedness of the sets A{x, ^) and A{x, —^). 

Lemma 5.11. If B2 is invertible and if at least one eigenvector of J is not an eigenvector of B, 
then for every x € R^ both sets A{x,^) and A{x, —^) are unbounded. 
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Figure 1: z belongs to ^(u),^a,u;) (the case r < 1). 

Proof. Take v € R"^ such that q{v) = —v x Bv — BJ^^{Jv x v) ^ O3; this is possible since at least 
one eigenvector of J is not an eigenvector of B. For every a G R the vector field x 1— > C{av)x + c{av) 
is in ^ and 

C(av)x + c(av) , , 
lim = q\v). 

This proves that, given a time T > 0, the solution of x = C{pLv)x + c(af) on the interval 
\—Ta~'^.,Ta~'^\ is as close as we want to the solution of i; = q{v) on [— T, T] (with the same ini- 
tial condition) provided that a is large enough. Now, the solution of x = g(v), x(0) = xq, leaves any 
fixed bounded set both in time T and in time — T, provided that T is large enough. □ 

We consider now, as in [15], the family {/ia,«, | A > 0, w G R"} of affine diffeomorphisms of R" 
defined by hx^w{x) = w + A(x — w). An easy computation shows that, for every X G 

{hx,^^X){x) = AiX){x - w) + XX{w) 

where h\^yj^X denotes the pushforward of the vector field X by /ia,u)- In particular 

lim(/iA,.^X)(x)=A(X)(x-«;). (26) 

A — >\j 

Denote by ^a,«) the image of the family ^ by h\^^^, that is, the transformation of ^ under the 
change of coordinates hx^w ■ 

Proof of Lemma \5.8l Let A{j^^) and r be as in the statement of Lemma 15.101 Consider the finite 
family A(^^) = { A(X)(x — w) \ X £ } and notice that each point of the sphere of center w and 
radius r is normally A(^^)-accessible from every point of the sphere of center w and radius 1. Let 
^\,w = {^A,u)^^ I ^ £ Thanks to formula (|26p we can assert that, if A is chosen sufficiently 
small, the sphere of center w and radius r is contained in the reachable set for ^a,w) from every point 
of the sphere of center w and radius 1 (see [211 Lemma 3.2]). Fix such a A > 0. 

Given w in R^ we claim that A{w,^) contains a neighborhood of w. Indeed, since the sets 
A{x, ^) and A{x, are unbounded, the same is true for the sets ^(x, ^\^w) and ^(x, — ^a,w))i 
for every x G R^. This implies in particular, because of the arc-connectedness of attainable sets, that 
there exists y G A{w., ^\^w) such that ||y — = 1. Let p < r and fix z G B{w., p). Since A{z, —^\w) 
is unbounded, then, again by arc-connectedness, there exists z' G A{^z^ —^\,w) such that tt;|| = r. 
Since from every point of the sphere of center w and radius 1, we can reach any point of the sphere 
of center w and radius r, the point z' is reachable from y by the family ,^\^w Finally, every point 
z G B{w,p) belongs to A{w,J^x^^) (see Figure 1). Since ^a,«) is the transformation of ^ by a 
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difFeomorphism preserving w, we deduce that A{w, contains a neighborhood of w. Therefore, 
every attainable set A{x, ^) is open. 

Since the control system defined by — =^ has the same form as that defined by ^ and verifies the 
hypotheses of both Lemmas 15. 101 and 15. IH the reasoning above proves also that every attainable set 
A{x^ —,^) is open. Take now y in the closure of A{x, The set A{y, being open, there exists 
z S A{y, — J^) n A{x, ^) which proves that y is reachable from x by the family ^ . Every attainable 
set A{x^ ^) is therefore both open and closed; thus it is equal to R'^. □ 

5.3 Generic controllability: physical interpretation 

The scope of this section is to provide a better physical insight of Theorem 15. li More precisely, 
we want to check that for m > 3 there exists an open and dense set in the (suitably defined) 
space of microscopic organisms modeled here such that the corresponding control system ©-([I]) is 
controllable. 

The space of organisms, denoted by Sm, will be identified with the set of pairs (0,^'), where 
$7 C is open, nonempty, connected, bounded, and of class C^, while \I' is the m-uple of functions in 
C^(9r2, R^) determining the action of the controls (see Section [3]). For simplicity we assume that the 
density 6 of the organism is constant, so that the inertia matrix J is determined (up to the positive 
multiplicative constant 5) by the shape $7. (For the general case see Remark l5.14l ) The topology on 

can be defined assigning the basis of open subsets defined by 

.yKn,e{^,-^)={{{ld + v){n),^) I V G W^^^{n,B^),^ G (C2(5fi,R3))™, 
||^^||3,oo < e, ll^* o (Id + v)- ^'||c2(Q) < e] 

for all £ G (0, 1/2) and all (17, ^) G (see [20] or, for a different approach, [8]). By extension Sq 
denotes the set of C^, open, nonempty, connected, bounded subsets of R^ endowed with the topology 
whose basis is given by all 

J^oA^) = {{ld + v){n) I V G T^3'~(f7,R3), ||t;||3,^ < e] , 

e G (0,1/2), O G So. 

Theorem 5.12. Assume that m > 3. There exists an open and dense set O in such that 

is controllable if (O, ^) belongs to O. Moreover, there exists an open and dense set V in Sq such 

that O n ({0} X (C2(90,R3))'^) is dense in {Q} x {C^ {0^,11^ for every n in V. 

Proof Denote by F{Q, ^) the triple {A, B, J) associated to an element (0, ^) of Sq through ^ 
and (jlOp . We want to show that each condition appearing in the statement of Propositions 15.21 and 
15.31 is satisfied by the elements of the image through F of an open and dense subset of Sq. 

Notice that the map ^' i— > i? is onto when considered from (C^(9r2, R'^))™ to ^exm fixed). 
Indeed, according to [191 Lemma 4.1], g^^\ g^'^\ g^^\ G^^\ G^'^\ G^^^ are linearly independent in 
L^(9r2,R3). Therefore, since the orthogonal to a smooth function in C'^{dQ,Ii^) with respect to 
the L^-product has infinite codimension, the relation (llOp defines a surjective map. In particular, 
the pre-image of an open and dense set in ^Qxm is open and dense in (C2(9Q, R^))*". 

Let us now take into account the dependence on We claim that the map from Sq to R^ 
that associates to a domain the spectrum of the corresponding inertia matrix -with eigenvalues 
repeated according to their multiplicity and with 6 normalized to one- is locally open at every point. 
In order to check it, fix 17 G Sq and a system of coordinates such that O3 is the center of mass of 
O and 61,62,63 its principal axes of inertia. Then the elements of the normalized inertia matrix 
jO = J/5 of n are 

tO _ / ± L ^i^jdx = if i j, , . 
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Apply the transformation Id + u to where v is the diagonal matrix diag(i;i, U2, fs), and denote 
by J^{v) the normalized inertia matrix of (Id + v){Q). Then 

70 ^,^ _ / ±(nLi(l + + ^0(1 + ^j) In XiXjdx = if i 7^ j, 

- I ±(n3^,(l + V,)) + v,fxl)dx if ^ = J. 

The spectrum of J^{v) is therefore given by the components of the vector a^{v) = (ib(n^_-^(l + 
Vk)) J^(Ylik^i{^ + ''JkYx'k)dx)^^i. A straightforward computation shows that the determinant of the 
Jacobian matrix of with respect to {vi,V2,vz), computed at vi = V2 = = 0, is different form 
zero. The map v 1— >■ a^(y) is therefore a submersion at t; = 0. 

In particular, the eigenvalues of J are simple for Q. in an open and dense subset of Sq, indepen- 
dently of ^. If 62 does not belong to any of the three planes generated by pairs of eigenvectors of J, 
then (jl7p is automatically verified. 

Similarly, one notices that all the assumptions appearing in the statements of Propositions 15.21 
and 15.31 can be represented in the form G{Vl, G O where O is an open and dense subset of a 
finite-dimensional vector space V , the map G : Sm ^ is continuous and for in an open and 
dense subset of Sq we have that G{Q., ^) belongs to O if belongs to an open and dense subset of 
•^dxm (possibly depending on 0). Therefore, G ^(O) is open and dense in 

Remark 5.13. As recalled above it follows from \19^ that, once is fixed, the linear map associating 
to the m-uple the matrix B through I is onto as a map from (L^((?r2, R'^))™ to >/#6xm- Thanks 
to this result, and to the remark that is controllable if m = 6 and B is invertible, the 

authors prove that for m > 6, for an open and dense subset of m-uples in C^(9r2,R^) (Q fixed) the 
corresponding system is controllable (see Theorem 1-1]). In the language adopted here the result 
of ]19^ says that, for every fixed $7 and for m > 6, controllability is a generic property with respect 
to 

Remark 5.14. In order to introduce a reasonable notion of genericity in the case where the density 
6 of the organism is not constant, it is necessary to include 6 in the definition ofT,m- Let be the 
set of triples {Q,,'^,5) where ($7,^) € and 6 S L°°{0,, {0,+oo)) and endow Tim with the topology 
whose basis is given by all 

jL,e{^, ^, 5) = { ((Id + v)(Vt),^, 7) I t; e I^3,oo(^^ r3), $ e {C\dVt, R3)r, 7 G L°°(r!, (0, +00)), 
||^^||3,oo < e, ll^* o (Id + v)- ^'||c2(n) < e, II7 ° (Id + v)- 5\\loo^^) < e} , 

e € (0,1/2), (r2,^,(5) S T,m- Following the same arguments of proof as above, one can prove that 
if m > 3 there exists an open and dense set O in such that is controllable if {Q,'^,6) 

belongs to O. The only difference consists in verifying that the map associating to (il, 5) the spectrum 
of J is locally open, and this can be done by taking the same perturbation (Id + diag(t;i, f2, U3))(r2) 
as above for Q and the perturbation 5 o (Id + diag(f 1, t;2, ^^3)) for 8. 

6 Swimming spherical microscopic organisms 

When the microscopic organism happens to be spherical, the equations presented in the previous 
sections have a very special form. As described in |19j . general results in hydrodynamics (see, e.g., 
|llj ) show that the matrix A corresponding to a spherical organism is diagonal and, more precisely, 
of the form 

A = ( -P'^"^ ° ^ 
V -p2\d ) ' 
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with 

P2 > Pi > 0. (28) 

Moreover, since the inertia matrix J is proportional to the identity (assuming that the mass distri- 
bution is homogenous inside 0), the nonhnear term E appearing in is given by 



E{z) 



03 



Notice that in the spherical case uj satisfies a well-defined control subsystem, ie, 

LO = —p2U) + B2U. (29) 

We are going to prove the following result. 

Theorem 6.1. Let Q be a ball and assume that the mass distribution of Vt is homogenous. Then 
the control system ^-^ controllable if and only if the rank of B2 is equal to three and Bi is 

different from Osxm- 

The easier part of the proof is to show that the controllability of ([S])^© implies that 

rank52 = 3, 5i / Osxm- (30) 

Indeed, the linear control system ([29]) is, as it is well known, controllable if and only if ranki?2 = 3. 
Moreover, if Bi = Osxm) then the space {O3} x Yi? is invariant for the dynamics of the control system 
([2]), which is therefore non-controllable. 

The converse implication will be proven in several steps. 

6.1 The case where Bi and B2 are hnearly independent 

In this section we study the case where the organism is spherical and 

rankSa = 3, Bi ^ span(S2). (31) 

Let us prove the following technical result. 

Lemma 6.2. // assumption h31\) holds true, then there exist ci, 02,03 G R^, an orthonormal basis 
{di, 1^2,^3} ofR?, and a m x 3 matrix T such that 

3 3 
Birv = '^ViCi, B2TV = '^Vidi, ci X di / O3. (32) 

i=l i=l 

Proof. Firstly, let us consider the case m = 3. Then B2 is invertible and for every orthonormal basis 
{di, (^2,^3} of R^, we have 

3 

B2^v = ^Vidi, 

i=l 

for r = i?2^^((ii|(i2|(i3), where {di\d2\d3) denotes the matrix whose columns are di, d2 and ^3. Since 

BiB2^ is not proportional to Id, then di can be chosen such that di x (iJiiJ^^di) 7^ O3. It suffices 

to complete di to any orthonormal basis {^1,^2,^3}, since BiTv = Ylf=i''^i'^i with q = SiS^^dj. 

(7) 

Let now m > 3. Denote by B^ the j-th column of B^. 

The second case we consider is when there exist 1 < Ji,j2>j3 ^ ^ such that B2^\ B2'^\ B2^^ 
are linearly independent and {b[^^^\b[^^^\b[^^^) span(5^^'^ l^^-'^^). Then, replacing v by 
T'v = {vj-i,Vj^,Vj^), we fulfill the hypotheses of the case m = 3 and we are done. 
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Finally, consider the case in which for every triple 1 < Ji,j2,j3 < "t- such that B2^\ -62"'^'', Bif'^^ 
are linearly independent, we have that {b[^^^\b[^''\b[^'^^) is proportional to {B!f^^\B^^''\B!f''^). Let 

{ju...Jk} = {j\l<j<m, Bi'^^03}. 

Then iB[^^^\- ■ ■ \b[^''^) is proportional to (-B^^'^l • • • \B^^''^). Assumption ([311) implies that there exists 
1 < i < ?^ such that b[''^ ^ O3 and B2^ = 03. It suffices then to replace v by T'v = {vj^ + 
Vj,Vj2, . . . ,Vjf,) in order to satisfy the hypotheses of the previous case. □ 

Notice now that for every Q G 50(3), the change of coordinates ^' = Q(,,u;' = QujX' = Cj-R' = 
RQ~^ preserves the dynamics of ©-(jj]). This is equivalent to say that, up to a change of coordinates, 
we can replace (i?i,i?2) by {QBi,QB2), for every choice of Q € 50(3). Let di,d2,ds be as in ([32]) 
and choose Q such that Qdi = Cj for f = 1, 2, 3 (where we recall that (ei)i=i^2,3 denotes the canonical 
basis in R^). The above transformations show that proving the controllability of the control system 
©-(III) under assumption (f3T]) is equivalent to proving it under the hypotheses that m = 3, B2 = Id, 
and that the first column of Bi is not proportional to ei. 

Denote by 61, b2, b-s the columns of Bi. Accordingly, the controlled vector fields are the three 
constant vector fields Xi = (6j,ej,06), i = 1,2,3. 

Lemma 6.3. // Ii31\) holds then system is Lie-bracket generating. 

Proof. Applying (0) in the case J = Id we obtain that for every vi,V2 € R'^ the constant vector 
field V = {vi,V2,0q) satisfies 

[[Xo,Xi],V] = {ei xvi-b,x 772,09). (33) 

In particular, taking i = 1 and V = Xi in ([33]) . we deduce that Z = {ei x 61,09) belongs to L. 
Applying again ([33]) to V = Z we get that, for i = 1,2,3, Zj = (e^ x (ei x 61), O9) belongs to L. 
Since Z, Zi, Z2, and Z3 span R^ x {O9}, it follows that L contains 

span(Xi,X2,X3) +R=^ x {O9} = R^ x {Oe}. 

As noticed in Remark 13. II this proves the lemma. □ 
The controllability of system ([2]) follows from Proposition 15.31 indeed. An = pild is symmetric 
and invertible and hypothesis (j3ip implies that Bi = B1B2 has a smaller set of eigenvectors than 
J = Id. Theorem 16. II is therefore proven, because of Proposition 14.1] under the stronger assumption 
that ([31]) holds true. 



6.2 The proportional case 

In order to complete the proof of Theorem 16.1] we have to prove that system (l2|)-(|3|) is controllable 
when 

rankS2 = 3, Bi = XB2, A / 0. (34) 

Using the same normalization argument as in the previous section, we can assume that m = 3 
and B2 = Id. Moreover, we can take A = 1 by performing the change of coordinates ^' = ^/\,uj' = 
oj, C = C/A, R' = R. Therefore, without loss of generality, Xi = {ei, 64,05) for i = 1,2, 3. 

Let us prove the Lie-bracket generating condition. We have 

[Xo,Xi] = {-piei,-p2ei,Rei,RS{ei)) 

and 

[[Xo,X,], [XQ,Xj]] = {06,2R{ej x e^),RS{ej x e^)). 
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A further computation yields 

[[[Xo,Xii [Xo,Xj]], [Xo,Xk]] = (06,3i?(efc x (e,- x e^)),RS{ek x (e, x e^))). 
Therefore, the Lie algebra L contains the following vector fields 
/eA f-piej\ / O3 \ 



O3 
VO3/ 



-P2ej 
\RS{ej)J 



O3 
2Rek 
\RSiek)J 



( O3 \ 

O3 

\RS{ei)/ 



for k,l = 1, 2, 3. As a consequence, L contains also 

/(p2-pi)eA / O3 \ 

(/9i - p2)ej 



O3 

Re-i 

\ RS{ei) J 



/03\ 
O3 

VO3/ 



/ O3 \ 

O3 
O3 

\RS{ei)J 

S0(3) (since 



for i,j,k,l = 1,2,3. These vector fields form a moving frame on the manifold 
Pi + P2)- 

In order to conclude the proof of Theorem 16. II we need to show that ([2]) is controllable. Applying 
Lemma 15.51 this turns out to be equivalent to the controllability of the following system, 



X = -piX - V X X + {p2 - Pl)v, 

with the control v taking values in R'^. Decomposing the control as v 



(35) 

wx + W with tf; G R and 
W orthogonal to x, equation (j35p rewrites as 

X = {-pi + {p2 - pi)w)x - W X X. 

Clearly, choosing = O3 allows to move as desired along span(x), while the choice w = pi / {p2 — pi) 
permits to attain the sphere of radius |[x|| centered at the origin. The controllability of ([2]) is thus 
proven and the proof of Theorem 16.11 completed. 



7 The case where the densities of the microscopic organism and of 
the fluid have the same order of magnitude 

As pointed out in [19\ Section 2] , the assumption that the densities of the microscopic organism and 
of the fluid have the same order of magnitude leads to stronger simplifications than in the general 
case studied in the previous sections. In particular, according to the remarks of page 6 the control 
system (2.18)-(2.27) presented in [19] reduces in this case to 

- Av + Vp = on F X (0, T), (36) 
divv = onFx(0,r), (37) 



viy, t) = Cit) + ivit) y<y + J2 Mt)My) 



i=l 



lim v{y, t) = 

\y\^oo 



an 



a{v,p)nds = 



an 



y X a{v,p)nds = 

C = RC 
R = RS{io) 



on F X (0,r), 
for {y,t) edUx (0,r), 
for t G (0,r), 

on (o,r), 

on (o,r), 

on (o,r), 
on (o,r). 



(38) 
(39) 

(40) 

(41) 

(42) 
(43) 
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Physically, v represents the field of velocities of the fluid in F (the coordinates are attached to the 
body), p is the pressure, ^ and u> are the linear and angular velocities of the organism and finally C 
and R give its position and orientation with respect to a fixed frame. Recall that . . . , are the 
functions characterizing the control actions and that a is the Cauchy stress defined in ([7]). 

Given a time-interval [0,T] and a control u € L°°([0, T], R™) there exists a unique solution 
{v, p, cj, C, R) of (IMD-dSD satisfying 

V G H^{[0, T],L'{F) n ^^''■(F) n D^^'^iF) n C°°(F)), 
supj^gj;'(l + < oo for almost every t S [0,T], 

p e H^{[0,T],L'-{F) n D'^^'^iF) nC^iF)), 
UH\[0,T],1l^), u;G//H[0,r],R3), C G CH[0, T], R^), ii G Ci([0, T], S0(3)), 

for s G (3,00], r G (3/2, 00] and 6 G (l,oo) (see [HIIH]). 

The finite-dimensional reduction of system (j36p - (j43p can be obtained following the procedure 
proposed in [19]. To this extent, as in Section [3l define the fundamental solutions (/i^*\p^*^) and 
{H^'''\ of the Stokes system and associate to them the matrices 0*^, T'^, k, /C. Let L be the nxm 
matrix 

91 02 
'^1 



L 



K 



Lemma 7.1. For almost every t G (0, T) we have 



Proof. According to |19^ Lemma 3.4] we have 



. y ds 

Ian 

G^'^ -yds = 



Lu{t). 



nds 



an 



Jan 

y X a{v,p)n ds 








where the equalities in the right-hand sides follow from ()40p and (|4ip . 
On the other hand, for almost every t G (0,T), 



• V ds 



an 



and, similarly 
Therefore, 



~ m 

/ 5^*) • (e + w X y + V -UjVj) = -{Q^C + + z^^^) 
Jan 

[ G^'^ ■vds = -{T^C + '^^^ + )Cu)i. 
Jan 



JC 



u 



01 92 



which proves the lemma. 

System ()36p -(|43 p therefore reduces to the control system 



□ 



C =RL^u, 
R =RS{L'^u), 



u G R" 



(44) 
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where L^,L^ G -^sxm are such that 



L 



L 



2 



By homogeneity, since the set of admissible controls is the whole R*", the controllability in 
arbitrary small time of system (|44|) is equivalent to the controllability of the control system whose 
admissible velocities are given by the family of vector fields = {±Xj | 1 < i < m} defined by 

X,{C,R) = {Rk,RS{ci)) (45) 

where the 6j's denote the columns of and the q's those of L^. Moreover, since J- is symmetrical, 
the controllability of this family is equivalent to the fact that the Lie algebra generated by denoted 
by L, is of rank 6 at every point of R'^ x S0(3) (see, eg, [U Corollary 5.11]). 



Theorem 7.2. The control system (j44p is controllable in arbitrary small time if and only if one of 
the following conditions is satisfied: 

1. rank(L2) > 2 and {L^'fL'^ + {L'^fL^ + 0„xm; 

2. rank(L2) = 2 and rank(L) > 3. 

Proof. First of all notice that the Lie bracket between two vector fields of the form Zi{(^,R) = 
{Rvi, R S{wi)), i = 1,2, is given by 

[Zi,Z2]{C,R) = {R{wi xv2-W2Xvi),RS{wi XW2)). (46) 

Let us prove the "if" part of the theorem. Suppose that we can find in L two vector fields Yi,l2 
of the form 

Y,{(:,R) = iRPi,RS{^i)), i = l,2, 

such that the vectors 71 and 72 are linearly independent and that the inner product (/3i,7i) is 
nonzero. We prove here below that, in this case, the Lie algebra L is of full rank at every point of 
R^ X S0(3). Then we show that both 1. and 2. guarantee that, without loss of generality, such Yi 
and Y2 can be found. 

Let Y3 = [Yi,Y2] whose expression can be obtained using (j46p . Define the vector field Wi as the 
bracket between I3 and Yi. Applying again (|46|) we obtain that 



WiiC,R) = [Y3,Yi]{C,R) 



i2(2(/3i, 71)72 - ((/?i,72) + (/32,7i))7i - (71, 72)/?! + ||7if /?2) 
i2 5(- (71, 72)71 + II71IP72) 



Let W2 be the vector field obtained adding (71,72) Yi — II71IPI2 to Wi, ie, 

W2{C,R) ■ 



i?(2(/3i, 71)72 - ((A, 72) + (A, 71)) 71; 
O3 

Since W2 is in L, then the vector field I4 = [W2,li] belongs to L as well and we have that 
YiiCR) = (i2(2(/3i, 71)71 X 72), O3). Consider now the vector fields Y5 = [Y4,Yi] and Yq = [Y4,Y2], 
both belonging to L, and whose expressions are Y5{(,R) = (i2(2(/3i, 71)71 x (71 X 72)),03) and 
Yq{C, R) = (i?(2 (/3i, 7i) 72 X (71 X 72)), O3). Since each matrix R in S0{3) is invertible, and be- 
cause (/?i,7i) ^ 0, the vector fields Yi,...,Yq span %,i?)R^ x 30(3) at every {C,R) G R^ x S0{3) 
if and only if the matrix 



A 



Pi P2 72 X /3i - 71 X (32 71 X 72 71 X (71 X 72) 72 X (71 x 72) 
71 72 72 X Ci O3 O3 O3 
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is full-rank. The non-degeneracy of A follows easily from the assumption that the vectors 71 and 72 
are linearly independent. 

We shall prove now that the hypotheses 1. and 2. imply the existence of two vector fields Yi and 
I2 as above. 

Let r = rank(L^). Notice that for every matrix T G GL{m) the reparameterization u — > Fu 
transforms the matrix {V-,L?) into {L^T,L?'T). The condition (L^) + (L^) / ^mxm is 
preserved by this transformation, since, {L^vfL'^T + {L'^TfL^T = {{L^f + {L^f L^)T . Thus, 
without loss of generality, we can assume that = (L^'^ lOsx (m-r)) where L^'^ G ^sxr is of rank r. 
Let us write as {L^'^ \ L^'"^) where L^'^ € ^sxr and L^'^ G -^3x{m.-r)- 

Suppose that 1. holds. If (L^'^) L^.i + (^2,1)^ ^1,1 ^ g^^xr, then we can find 1 < Ji,j2 < t such 
that {bj^,Cj^) ^ 0, or {bj^,Cj.,) / 0, or {hj^,Cj.,) + {hj^,Cj^) 7^ 0. If {bj^,Cj^) = {hj^,Cj^) = 0, then 
the inner product {hj^ + bj^^Cj^ + Cjj) = (bjnCj^) + {bj^^Cj-^) is nonzero, and thus the vector fields 
Yi = Xj^ + and Y2 = Xj^ fulfill the required conditions. Otherwise, {L'^'^) L^'^ ^ Orx{m-r) 
and we can fix 1 < ji,j2 < r and js > r such that ji / j2, {bj^,Cj^) = and {bj.^,Cj-^) 7^ 0. Then 
Yi = Xj^ + Xj^ and Y2 = -'^j^ fulfill the required conditions. 

Assume now that 2. holds and 1. does not. Then ci and C2 are linearly independent and 

(61, ci) = (62, C2) = (61, C2) + (62, ci) = 0. 

Moreover, there exists j > r = 2 such that bj 7^ O3 and {bj, ci) = {bj, C2) = 0. Therefore, there exists 
a real number a 7^ such that 63 = a(ci x C2). Then, according to (j46p . 

[XuX,]{C,R) = {R{aci X (ci X C2)),03). 

As (c2,ci X (ci X C2)) 7^ the vector fields Yi = Xi and 5^2 = ^^2 + fulfill the required 

conditions. 

Let us prove now the "only if" part of the statement. First notice that if r < 1 then there exists 
a vector xq € R'^ \ {O3} such that L'^u x xq = for every control u. Multiplying equation (j44p by 
xo, we get that Rxq = and so R{t)xo = R{0)xo for every time t. Therefore, system (j44l) is not 
controllable. 

Assume now that r = rank(L) = 2 and (L^) + (L^) = Omxm- Let Y = [Xi,X2], whose 
expression, according to (|46p . is given by 

Y{C,R) = {R{ci X 62 -C2 X bi),RS{c^ x C2)). 

Since Cj) = for i = 1, 2 and (61, C2) + (62, ci) = 0, we obtain that 

\Y,Xi] = -(ci,C2)Xi + ||ci||2X2 
[y.X2] = -||c2fXi + (ci,C2)X2. 

Therefore, the Lie algebra L, which is generated by Xi and X2, is equal to the linear space of vector 
fields spanned by Xi, X2 and Y and cannot be of full rank. 

Let now r = 3 and (L^) + (L^) L} = Omxm- The condition [L"^'"^)^ L^''^ = 03x(m-3) implies 
that the columns of L^'^ are orthogonal to all the elements of a basis of R^. Therefore, L^'^ = 
03x(m-3)- As for the columns of L^'^, we easily obtain from {L^'^) Lp"^ + [L"^'"^) L^'^ = O3X3 that 

61 = - ai2 C2 - ai3 C3, 

62 =ai2Ci - 023 C3, 

63 =Qi3 ci +023 C2, 
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for some 012,013,023 ^ Without loss of generality we can assume that (01,02,03) is a positively 
oriented orthonormal basis and an easy computation gives 



which proves that L is equal to the linear space of vector fields spanned by Xi, X2 and X^ and 
cannot be of full rank. □ 
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